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Abstract 

Gauge/gravity duality is a useful and efficient tool for addressing and studying questions 
related to strongly interacting systems described by a gauge theory. In this manuscript we 
will review a number of interesting phenomena that occur in such systems when a background 
magnetic field is turned on. Specifically, we will discuss holographic models for systems 
that include matter fields in the fundamental representation of the gauge group, which are 
incorporated by adding probe branes into the gravitational background dual to the gauge the- 
ory. We include three models in this review: the D3-D7 and D4-D8 models, that describe 
four-dimensional systems, and the D3-D7' model, that describes three-dimensional fermions 
interacting with a four-dimensional gauge field. 



*To appear in Lect. Notes Phys. "Strongly interacting matter in magnetic fields" (Springer), edited by D. Kharzeev, 
K. Landsteiner, A. Schmitt, H.-U. Yee. 
tbergman@physics.technion.ac.il 
'tjke@mppmu.mpg.de 
§ giladl @research.haifa. ac.il 



1 



1 Introduction 



The behavior of strongly interacting matter subject to background magnetic fields is an interesting 
and physically relevant problem in many different scenarios, ranging from the effective 2d electron 
gas in graphene, to magnetars, which are neutron stars with a strong magnetic field. Magnetic 
fields give rise to a rich array of phenomena. Some examples in QCD are the magnetic catalysis of 
chiral symmetry breaking 03 El El, anomaly-driven phases of baryonic matter H, and the chiral 
magnetic effect |0. It has also been suggested that magnetic fields induce p -meson condensation 
and superconductivity in the QCD vacuum [|6lH- There are also many interesting examples in 
condensed matter physics, most notably the fractional quantum Hall effect |H|9l- 

Gauge/gravity duality, also known as holographic duality, has emerged in recent years as a par- 
ticularly useful approach to strong-coupling dynamics. Although it does not seem to be directly 
applicable to physical systems, this approach can be used to study theoretical systems that exhibit 
the same type of phenomena, and that capture some of the relevant physics. The techniques of 
holographic duality are especially efficient in addressing questions associated with finite tempera- 
ture and density, background fields and transport properties, that are difficult to study using other 
non-perturbative methods. This approach can also lead to new ideas for constructing effective 
theories of the physical phenomena one is interested in. 

Holographic models are divided into two main types, commonly referred to as top-down mod- 
els and bottom-up models. In top-down models the bulk gravitational description of the system 
corresponds to a consistent solution of a well-defined quantum gravity theory, either in the con- 
text of the full string theory or in terms of the low-energy effective supergravity theory. This then 
defines some particular strong-coupling boundary dynamics. In bottom-up models, on the other 
hand, one builds into the description what one needs in order to produce the desired boundary 
dynamics. Each approach has advantages and disadvantages. Top-down models are more firmly 
grounded than bottom-up models, however they are more restrictive in terms of the variety and 
scope of phenomena they can exhibit. 

Probe-brane models are a class of top-down holographic models, in which matter fields trans- 
forming in the fundamental representation of a gauge group are incorporated by embedding "fla- 
vor" D-branes in the gravitational background dual to the gauge theory IfTOl . These branes are 
treated as probes, in the sense that we neglect their backreaction on the background. (This cor- 
responds to the 'quenched' approximation in the dual gauge theory, where matter loops are ne- 
glected in computing gluon amplitudes.) The matter fields are manifest in this construction: they 
correspond to the open strings between the flavor branes and the "color" branes that make up the 
background. In particular, one can easily design probe brane models in which the light matter 
degrees of freedom are purely fermionic, which is obviously a desirable feature for many physical 
systems, including QCD and condensed matter electron systems. 

The matter dynamics is determined by the properties of the probe brane embedding. In particu- 
lar, the fluctuations of the probe brane worldvolume fields correspond to gauge-invariant composite 
operators that describe the mesonic states of the matter system. There are generically two types of 
embeddings at finite temperature: "BH embeddings", in which the brane extends to the horizon of 
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the background, and "MN embeddings", in which the brane terminates outside the horizon. The 
two embedding types describe different phases of the matter in the dual gauge theory. For example, 
in the MN phase the mesons are stable since they are associated with real eigenfrequencies of the 
probe brane fluctuations. In the BH phase, on the other hand, some of the energy of the fluctua- 
tions is dissipated into the black hole, leading to complex eigenfrequencies and damping. In this 
case the mesons have a finite lifetime. MN embeddings are favored at low temperature, and as the 
temperature is increased one generically observes a first order phase transition to a BH embedding. 

The most extensively studied probe-brane models are the D3-D7 model IfTOl , in which D7- 
branes are added to the D3-brane background, and the D4-D8 (or Sakai-Sugimoto) model |[TT|. in 
which D8-branes and anti-D8-branes are added to the background of D4-branes compactified on 
a circle. Both models describe a strongly-coupled gauge theory in four dimensions with funda- 
mental matter degrees of freedom, and both exhibit a number of phenomena similar to QCD. More 
recently, a different D3-D7 system, more closely related to the D4-D8 system, has been used as a 
model of strongly-interacting fermionic matter in three spacetime dimensions [fT2l[T3l[T4l[T5l . We 
will refer to this as the D3-D7' model. This model exhibits several interesting phenomena that are 
familiar in planar condensed matter systems. 

Probe-brane models are especially well-designed to study the properties of the dual matter 
systems at non-zero density and in background electromagnetic fields. Both are implemented 
by turning on specific components of the probe-brane worldvolume gauge field, and solving the 
resulting coupled differential equations for the embedding and the gauge fields. Here one observes 
another basic difference between the two types of embeddings in terms of their response to a 
background electric field. MN embeddings correspond to electrical insulators with a mass-gap to 
charged excitations, and BH embeddings describe gapless conductors. 

Probe-brane models also exhibit a number of interesting phenomena in a background magnetic 
field, which are qualitatively similar to the phenomena listed in the beginning. In this paper we will 
review how each of the three models mentioned above respond at non-zero density to a background 
magnetic field in various situations. In particular, we will encounter the magnetic catalysis effect in 
both the D3-D7 and D4-D8 models. In the D3-D7 model we will also demonstrate the formation of 
a superfluid state. In the D4-D8 model we will describe anomaly-generated currents and baryonic 
states, as well as a metamagnetic-like transition. In the D3-D7' model we will see both quantum 
and anomalous Hall effects, as well as how the magnetic field influences the instability to the 
formation of stripes, and the zero-sound mode. 

This paper is divided into three main sections, reviewing each of the probe-brane models in 
turn. For completeness let us mention that magnetic fields also play an important role in the related 
D3-D5 system, where the probe D5-brane corresponds to additional (2+ 1) -dimensional degrees 
of freedom in the dual gauge theory. In this model, the magnetic field leads to a phase transition of 
Berezinskii-Kosterlitz-Thouless (BKT) type [fT6l ITVTl . For brevity we do not discuss this model in 
this review. 
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2 The D3-D7 model 



2.1 Brane construction 

The starting point for this model is the usual configuration of the AdS/CFT correspondence |fT8ll 
which involves a stack of N D3 branes. This has an open string interpretation in which the low- 
energy degrees of freedom are described by U (N) JV = 4 Super- Yang-Mills theory. On the other 
hand, in the closed string interpretation of N D3 branes, the low-energy near-horizon limit gives 
rise to the space AdS$ x 5 5 . Identifying the two pictures leads to the AdS/CFT correspondence. 

Let us now add Nf probe D7-branes to this configuration, as first done in |[T0l and reviewed 
in detail in [fT9l . Within (9+l)-dimensional flat space, the D3-branes are extended along the 0123 
directions, whereas the D7-branes are extended along the 01234567 directions. This configuration 
preserves 1 /4 of the total amount of supersymmetry in type IIB string theory (corresponding to 8 
real supercharges) and has an 50(4) x 50(2) isometry in the directions transverse to the D3-branes. 
The 50(4) rotates x 4 ,x 5 ,jc 6 ,jc 7 , while the SO(2) group acts on jc 8 ,jc 9 . Separating the D3-branes 
from the D7-branes in the (8,9) directions by a distance I explicitly breaks the 50(2) group. These 
geometrical symmetries are also present in the dual field theory: The dual field theory is an jY = 2 
supersymmetric (3+l)-dimensional theory in which the degrees of freedom of jY = 4 super Yang- 
Mills theory are coupled to Nf hypermultiplets of flavor fields with fermions and scalars (v?',? n )> 
i = 1 , 2, n = 1 , 2, which transform in the fundamental representation of the gauge group. Separating 
the D7-branes from the D3-branes corresponds to giving a mass to the hypermultiplets. 

For massless flavor fields, the Lagrangian is classically invariant under conformal transforma- 
tions 50(4,2) F] Moreover, the theory is invariant under the ^-symmetries SU(2)r and U(1)r as 
well as under the global SU (2)$, which rotates the scalars in the adjoint hypermultiplet. Note that 
the mass term in the Lagrangian breaks the U(1)r symmetry explicitly. If all Nf flavor fields have 
the same mass m, the field theory is invariant under a global U (Nf) flavor group. The baryonic 
U(1)b symmetry is a subgroup of the U (Nf) flavor group. These symmetries of the field theory 
side may be identified with symmetries of the D3-D7 brane intersection and hence also with the 
dual gravity description. 

For this field theory, gauge invariant composite operators may now be constructed which trans- 
form in suitable representations of the SU (2) x SU (2) x U(l) symmetry group isomorphic to the 
geometrical 50(4) x 50(2). These operators are expected to be dual to the fluctuations of the D7- 
brane which transform in the same representation, as worked out in detail in ||20l . An example of 
a meson operator is given by 

J( A = Xjfi <7 A lj y/j + q m X A q m , (/, m = 1, 2) , ( 1 ) 

with X A the vector (X s , X 9 ) of adjoint scalars associated with the (8,9) directions, and o A = 
(a 1 , a 2 ) a doublet of Pauli matrices. Thus ([T]) has charge +2 under U(1)r. It is a singlet under 

'However note that the scale-invariance is broken at the quantum level since the beta function is proportional to 
Nf/N c and therefore non-vanishing. In the limit N c —> °° with Nf being fixed, the beta function is approximately zero, 
i.e. we may treat the theory as being scale invariant also at the quantum level. 
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both SU(2)<$ and SU(2)r. The conformal dimension is A = 3. This operator may be viewed as 
a supersymmetric generalization of a mesonic operator in QCD, with the index A labelling two 
scalar mesons. 

The standard AdS/CFT duality relates the jY = 4 Super Yang-Mills degrees of freedom to 
supergravity on AdSs x S 5 . In addition, there are new degrees of freedom associated to the D7- 
brane worldvolume fields originating from the open strings on the D7-brane. The additional duality 
maps these to the mesonic operators in the field theory. This is an open-open string duality, as 
opposed to the standard AdS/CFT correspondence, which is an open-closed string duality. The 
dynamics of the D7-brane is described by the Dirac-Born-Infeld (DBI) action 

Sdi = -—J d 8 ^-det (Gab+B^+lna'Fa^ , (2) 

where /I7 = [(27r) 7 G5 /4 ] _1 . G and B^> are the induced metric and two-form field on the probe 
brane worldvolume, and F a h is the worldvolume field strength. The D7-brane action also contains 
a fermionic term S D1 . In addition there may also be contributions of Wess-Zumino form. An 
example for this will be discussed below. 

Let us write the AdS$ x S 5 metric in the form 

r 2 R 2 

ds 2 = —^r\ijdx l dx j + -^[dp 2 + p 2 dQ.l + dxl + dxj) , (3) 
R l r L 

with i, j = 0, 1, 2, 3, p 2 = x\ + ... + x 2 , r 2 = p 2 + x| + x\ and R the AdS radius. Since the D7-brane 
is transverse to x$, xg in flat space, we see that it extends along AdSs and wraps an S 3 inside S 5 in 
the near-horizon background. The action for a static D7-brane embedding, for which F a t, may be 
consistently set to zero on its world-volume, is given from up to angular factors by 

Sdi = -— j d^p 3 ^/l+x 2 +x 2 , (4) 

where a dot indicates a p derivative (e.g. is = d p x%). The ground state configuration of the D7- 
brane then corresponds to the solution of the equation of motion 

d 

dp 

where x denotes either jcs or xg. Clearly the action is minimized by x%,xg being any arbitrary 
constant. Therefore the embedded D7-brane is flat. According to string theory, the choice of the 
position in the x%,x<) plane corresponds to choosing the quark mass in the gauge theory action. 
The fact that ^8,^9 are constant at all values of the radial coordinate p, which corresponds to the 
holographic renormalization scale, may be interpreted as non-renormalization of the mass in the 
dual field theory. 

In general, the equations of motion have asymptotic (p — > 00) solutions of the form 

X=l+Ay+..., (6) 
p 2 



p dx 

1 ~\~ Xq ~\~ Xq ^ 



0, 



(5) 
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where / is related to the quark mass m by 



In agreement with the standard AdS/CFT result about the asymptotic behaviour of supergravity 
fields near the boundary, the parameter c must correspond to the vev of an operator with the same 
symmetries as the mass and of dimension three, since p carries energy dimension, c is therefore a 
measure of the quark condensate i/ai/a. c is obtained from /dm which in addition to the fermion 
bilinear also includes scalar squark terms. We may consistently assume that the squarks have zero 
vev. Moreover, supersymmetry requires that a vev for c must be absent since c is an F-term of a 
chiral superfield: ty/ty/ is the F-term of QQ. Supersymmetry is broken if c = (fyy) 7^ 0. This is 
reflected also in the supergravity solution: The solutions to the supergravity equations of motion 
with c non-zero are not regular in AdS space and are therefore excluded. 

We therefore consider the regular supersymmetric embeddings of the D7-brane for which the 
quark mass m may be non-zero, but the condensate c vanishes. For massive embeddings, the D7- 
brane is separated from the stack of D3-branes in either the x% or xg directions, where the indices 
refer to the coordinates given in ([3]). In this case the radius of the S 3 becomes a function of the 
radial coordinate r in AdS 5. At a radial distance from the deep interior of the AdS space given by 
the hypermultiplet mass, the radius of the S 3 shrinks to zero. From a five-dimensional AdS point 
of view, this gives a minimal value for the radial coordinate r beyond which the D7-brane cannot 
extend further. This is in agreement with the induced metric on the D7-brane world- volume, which 
is given by 

p2_|_^2 ^2 R 2 p 2 

ds 2 = _ riijdx'dxj + -= — ^dp 2 + = - in dQ. 2 , (8) 
R 2 J p 2 + l 2 p 2 + l 2 

dQ. 2 = d\]/ 2 + cos 2 yd [5 2 + sin 2 y/dy 2 , (9) 

where p 2 = r 2 — I 2 and ^3 are spherical coordinates in the 4567-space. For p — > °°, this is the 
metric of AdSs x S 3 . When p = (i.e. r 2 = I 2 ), the radius of the S 3 shrinks to zero. 



2.2 Finite temperature 

The finite temperature system is realized holographically by placing the D7-brane in an AdS- 
Schwarzschild black hole background with metric given by 




where r 2 = p 2 + L 2 . The temperature is given by T = r^j (nR 2 ). 
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In this metric we have introduced polar coordinates (L, 0) in the (x^^xg) plane and consider 
solutions for D7-brane embeddings with L = L(p), = const. The asymptotic near-boundary 
behaviour of these brane embeddings is given by 

L(p)=Z + -^ + ..., (11) 

where m = l/(2na') is the bare quark mass and c is proportional to the condensate (\jf\jf). At 
finite temperature, supersymmetry is broken and brane embeddings with c non-zero are possible, 
in contrast to the supersymmetric case discussed above. 

Depending on the quark mass, there are two different types of embeddings: Those which reach 
the black hole, and those which do not since the S 3 they wrap shrinks to zero outside the black 
hole horizon. The first type of branes is referred to as 'black hole' (BH) embeddings, while the 
second type is referred to as 'Minkowski' (MN) embeddings. In the BH case, fluctuations of the 
probe brane have complex eigenfrequencies or quasi-normal modes, which means that the mesons 
associated with these fluctuations decay. In the MN phase, the mesons are stable. The phase 
transition between the two types of embeddings is first order. 

At finite temperature, the solution with m = also has c = 0. However, in gravity backgrounds 
corresponding to confining field theories, brane embeddings with m = 0, c ^ are possible. These 
realize spontaneous chiral symmetry breaking [1211 . 



2.3 Magnetic catalysis 

Let us now consider, as was first done in (|22|. a magnetic field induced by a pure gauge B-field in 
the worldvolume direction of the D3-branes, 

5 (2) =Bdx 2 Adx 3 , (12) 

which satisfies dB^ = 0. This contributes to the DBI action ([2]). Since B^ and Ithx'F enter (|2]) in 
the same way, we may trade B^ for a gauge field on the probe brane via F = —B^ /2na', which 
justifies interpreting B of ( fT2] ) as a magnetic field. 

In addition, there is a non-trivial Wess-Zumino contribution to the action, which at first order 
in a' is of the form 

Swz = 2na' ii 7 Jf A C (6) . (13) 

In the presence of the B-field, this leads to an additional non-trivial contribution to the action, as 
explained in detail in 11221 . This gives rise to a non-trivial which breaks supersymmetry on the 
worldvolume of the D7 -brane. 



Since supersymmetry is broken, the D7-brane now has a profile which depends on p as in ( |TT[ ), 
even at zero temperature. The Lagrangian corresponding to ([2]) takes the form 



if = -^Vsinycos l/Vl +L /2 Jl + t~2_ A + L iy • ( 14 ) 
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For m = 0, the brane embedding solution obtained from this Lagrangian has non-zero c oc ({fry/) in 
( fTTj ). The magnetic field therefore induces spontaneous chiral symmetry breaking, a phenomenon 
known as magnetic catalysis [Q~l[2l[3]|. For large m, the condensate may be calculated analytically 
[|22l and is found to be 

(W)~-c = ~B 2 . (15) 

For small m, c has to be evaluated numerically. The result is shown in figure [T] By evaluating the 
free energy, it has been shown that when there is more than one solution, the one with the larger 
condensate is preferred. 




At small values of the magnetic field, it is possible to analytically evaluate the shift of the 
meson masses due to its presence. For fluctuations of the embedding scalar in particular, a Zee- 
man splitting is observed fl22). While in the absence of a magnetic field, in the supersymmetric 
case described here, the scalar meson mass obtained from the fluctuations is Mq(h) = 2m/\fX ■ 
\/(n+ l)(« + 2) [|20l . for non-zero magnetic field there is a mass splitting 

M±=M ±4f-- (16) 

A review of magnetic catalysis in probe D7-brane systems is given in ll24ll . Magnetic catalysis 
is also found in systems involving Nf D7 branes where the backreaction of the metric on the 
background geometry is taken into account [|25ll26ll27l . Out-of-equilibrium dynamics associated 
with the phase transition induced by magnetic catalysis has been investigated in ||28l . 

In the finite temperature case, there is a competition between two mechanisms: The black hole 
attracts the D7-brane, while it is repelled at small radii by the magnetic field. This implies a phase 
transition between a phase where the D7-brane reaches the black hole and one where it does not. 
This is shown in figure [2] for different values of the magnetic field, where the dimensionless ratio 
B/T 2 is used. A detailed discussion of the normalization is found in ||23l . 

There is a critical value for B/T 2 above which the probe brane is repelled from the black hole 
for all values of the bare quark mass: In this case, chiral symmetry is spontaneously broken and 
the mesons are stable. The phase diagram is shown in figure [3] 
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Figure 2: Increasing values of B/T 2 for fixed T show the repulsive nature of the magnetic field. 
We see that for large enough B/T 2 , the melted phase is never reached, and the chiral symmetry is 
spontaneously broken. Figure reproduced from lf23l . 



B 

T 2 

15 - 




Figure 3: Phase diagram for the D3/D7 system in the (B/T 2 ,T /m) plane. Figure from [|23l . 

More involved phase diagrams are obtained if a U( 1) chemical potential and density are turned 
on in addition to the magnetic field by considering a non-trivial profile for the U( 1) gauge field on 
the D7-brane ll29l[30ll3Ti . An example of a phase diagram is shown in figure |4j 

2.4 Superfluid 

At finite isospin density, the D3-D7 model realizes a holographic superfluid [|32l[33l . Finite isospin 
density is obtained by considering two coincident D7-branes, and using an ansatz for solving the 
equations of motion, which involves a non-trivial profile for the the temporal component of the 
SU (2) worldvolume gauge field, with asymptotic behavior 

A 3 t (p)oc^ + ^. (17) 
P 2 
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T 




Figure 4: Phase transitions at constant B field in the (T,ji) plane for the D3/D7 model. Figure 
reproduced from OTTl by kind permission of the authors, c, d, J refer to the condensate, density 
and electric current, respectively. 



/l 3 breaks the SU (2) symmetry explicitly to a residual £7(1)3. In the presence of this background, 
the energetically favored solution also involves a non-trivial spatial component of the worldvolume 
gauge field, 

P 

Here the leading contribution is absent in the asymptotic behavior, so the £7(1)3 symmetry is 
spontaneously broken. A l x is dual to a condensate of the form 

d 1 x °c{yo l y x }ir+$o i d x Q), (19) 

which is the supersymmetric equivalent of the p meson. The calculation of the frequency-dependent 
conductivity <j(ft>) for this solution shows that it describes a superfluid: o(co) displays a gap. For 
the Sakai-Sugimoto model discussed below, a similar condensation mechanism has been found in 
[|34l and superfluidity has been discussed in 11351 . 

As discussed in 1136113711 . a similar condensation process also happens when the profile ( fTTj ) for 
the temporal component of the SU (2) gauge field is replaced by a non-trivial profile for a spatial 
component of the form 

A\ = By, (20) 

which corresponds to a background magnetic field. In this case a similar condensation as above 
takes place. This has been demonstrated by analyzing the fluctuations about the magnetic field 
background ll37ll : The quasi-normal modes of particular fluctuations cross into the upper half of 
the complex frequency plane above a critical value of the magnetic field, indicating an instability. 
This is shown in figure [5} Finally let us note that the Hall conductivity has been calculated for 
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the D3/D7 model in 113811 . Unlike the isospin case, the ground state involving the p condensate 
is spatially modulated for the magnetic field background, leading to an Abrikosov lattice ||39ll . A 
similar p meson condensation mechanism in a background magnetic field has been found in the 
context of field theory in [|6l|71|40]|, based on similar earlier results in electroweak theory [j4TT| . For 
the Sakai-Sugimoto model which we discuss below, a similar mechanism has been discussed in 

num. 



Im ro 



Re ro 




Figure 5: Quasi-normal modes cross into the upper half plane above a critical magnetic field, 
signalling an instability. Figure reproduced from |[37ll . 



3 The D4-D8 (Sakai-Sugimoto) model 

3.1 Basics 

N c D4-branes on M 1 ' 3 x S 1 with anti-periodic boundary conditions for fermions provide a holo- 
graphic model for the low energy behavior of 4d SU (N c ) Yang-Mills theory with g\ M = ATtg s \fa/ /R4 
(|44|. The near-horizon background at zero temperature is given by (we work with dimensionless 
coordinates rescaled by R) 

ds 2 con = w2 (-dxl + dx 2 + f(u)dxj) +u~i \^—^+u 2 dQ?A 

e * = gsU 3/4 , F 4 = 3n(a , ) 3/2 N c da 4 , (21) 

where f(u) = 1 - (u 3 KK /u 3 ), u KK = 4R 2 /(9R 2 4 ) and R = (7tg s N c ) 1/3 Va*. The IR "wall" at u = u KK 
implies that the dual gauge theory is confining. At nonzero temperature there is another possible 
background with a metric 

3 3 / du 2 \ 

ds 2 dec = ui (-f(u)dx 2 ) + dx 2 + dxj)+u-2 \ -^+u 2 d£l 2 A , (22) 
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where f(u) = 1 — (uj/u 3 ) and uj = (4n/3) 2 R 2 T 2 . This background becomes the dominant one 
when T > 1 / (litRa) . The presence of a horizon at u = uj in this background indicates that the 
gauge theory undergoes a (first order) deconfinement transition at this temperature. 

Quarks are added to the model by including D8-branes and anti-D8-branes that are localized on 
the circle |fTT|. With Nf D8-branes at one point and Nf anti-D8-branes at another point, the model 
has Nf flavors of massless right-handed and left-handed fermions, and a U (N/)r x U (Nf)t chiral 
symmetry. The 8-branes are treated as probes in the near horizon background of the D4-branes. 
The flavor dynamics is thus encoded in the 5d effective worldvolume theory of the D8-branes, 
which includes a DBI term and a CS term (in Lorentzian signature^ 

Sdbi = J d 4 xdu u 1 / 4 \/-det(gMN + Imn) (23) 

Scs = ~Jd 4 xdue MNp Q R a M fNpf Q R. (24) 

The dimensionless worldvolume gauge field cim and field strength /mn are defined as % = 
illta' /R)Am and f MN = 2%a!F MN , and the overall normalization is given by JV = ^ A R 9 /g s = 
(l/3)N c R 6 (27l)~ 5 ((x')~ 3 . The anti-D8-brane has a similar action in terms of its worldvolume 
gauge field clm- The DBI term is identical to that of the D8-brane, and the CS term has the op- 
posite sign. We define the vector combination as a v M = j(aM + dM), and the axial combination as 

A 1 1 - \ 
a M = 2 \ aM ~ a M)- 



In the low-temperature confining background (21 ) the D8-brane and anti-D8-brane connect at 
u = uq> ukk into a smooth U-shaped configuration (Fig. [6^), reflecting the spontaneous breaking 
of the U(1)rXU(1)l chiral symmetry to the diagonal U(l)y. The embedding is determined by 
the DBI action (setting f MN = 0) 



$™ 1 = -jY J d 4 xduu 4 
which implies an asymptotic behavior 



f{u){x' A {u)f + 1 



u 3 f(u) 



(25) 



X4(u)^-~- u9/2 , (26) 

where L is the asymptotic brane-antibrane separation. 

The normalizable fluctuations of the D8-brane worldvolume fields in this embedding corre- 
spond to the (low spin) mesons of the model. Their mass scale is set by uq, which we can think of 
as the mass of a "constituent quark" described by an open open string from uq to ukk- There is one 
massless pseudoscalar field <p, precisely as one expects from the broken chiral symmetry. This is 



2 For simplicity, we will consider the single flavor case with one D8-brane and one anti-D8-brane. This does not 
affect any of the results qualitatively. 
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related to the 77 ' meson in QCD. It appears (in a gauge with a u = 0) as the zero mode a^: 

a^(x^,u) = - d^i (p (V 1 ) y/b ( u) + higher modes, (27) 

where 




y/b(w) = - arctan . / 1 . (28) 



U KK 



Baryons are described by D4-branes wrapped on S 4 inside the D8-brane. Their charge comes 
from the N c strings which must be attached to the wrapped D4-brane to cancel a tadpole due to the 
background RR field. These strings end on the D8-brane, giving iV c units of charge]^] 



In the high-temperature deconfining background (22) the D8-branes and anti-D8-branes can be 
either connected (Fig. [6Jd) or disconnected, with x${u) = L/2 (Fig. |6j:), the latter corresponding to 
the restoration of the chiral symmetry [|46ll . The DBI action in the high-temperature deconfining 
background is very similar: 



S d D % = -jK J d 4 xduu 4 



f(u){x' A {u)) 2 + 1 



u 3 



1 



(29) 



Consequently the properties of the U embedding in this background are qualitatively similar to 
those of the embedding in the confining background, for example in terms of the spectrum of 
mesons. In the disconnected embedding there are no normalizable fluctuations corresponding to 
mesons, as one expects in a chiral- symmetric phase. Comparing the (Euclidean) actions of the two 
embeddings shows that the disconnected one becomes dominant when T > 0. 154/ L. In particular, 
for small L(L< 0.97i?4) the gauge theory has an intermediate phase of deconfinement with broken 
chiral symmetry. 



3.2 Finite density and background fields 

The D8-brane worldvolume vector and axial gauge fields are dual to conserved vector and axial 
currents in the gauge theory, and therefore^] 

on— shell /ir>\ 

^V4da]r{u^°°) 

3 Note that, although the boundary value is non-zero, this is a normalizable mode since the field strength is nor- 
malizable. Ordinarily, boundary values of bulk fields correspond to parameters in the boundary theory. But in this case 
there is a possible ambiguity, since the boundary value of can also describe a non-trivial gradient of the pseudoscalar 
field. 

4 For Nf > 1 the baryons correspond to instantons in the non-abelian D8-brane theory lfTTl|45l . This reproduces 
the known description of baryons as Skyrmions in the chiral Lagrangian. In this description the baryon charge comes 
from the CS term coupling the U(l)y field to the instanton density in the SU(N/)v part. 

5 The axial symmetry is broken by an anomaly. However this is a subleading effect at large N c which we can 
neglect. In particular, we will assume that the one-flavor pseudoscalar rj' is massless. For a discussion of the U(1)a 
anomaly and the 77' mass in the context of the Sakai-Sugimoto model see lfTTll47l . 
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The chemical potentials are defined bjQ 

/ly = <Jq (u — > oo) and jU^ = Oq(w — )■ °o) . (31) 

In our conventions quarks carry one unit of vector charge and baryons carry N c units. Neverthe- 
less we will refer to the vector current as the "baryon number current". We are also interested 
in studying the effects of background "electromagnetic" fields that couple to this current, which 
correspond to turning on spacetime dependent boundary values of the worldvolume gauge field, in 
particular 

«i = foi(u -»■ °°) , h = £ijkfjk{u -»■ oo) . (32) 

In some situations one may be required to add boundary terms to the action. These are es- 
pecially relevant if there is a CS term in the bulk. In deriving the equations of motion from the 
variational principle one usually assumes that the surface terms vanish. However in some instances 
one has to be more careful. The surface terms (in the a u = gauge) are given in general by 



8S 



. a x -j-j- oa v 

on—shell J OCl v 



+ I a xdu^r-p: r oa v 

d(d^a v ) Xfl - 



(33) 



In holography the boundary values of the fields at u — > °° are fixed, so da^ (w — > °°) = 0. However 
8a^{u m i n ) and 8a^(x jl — > ±°°) need not vanish. Therefore a surface term may be non-trivial if 
the fields extend to these boundaries. In order to have a well-defined variational principle we must 
therefore add boundary terms S^{u m i„) +S^(x^ — > ±°°), whose variation cancels the surface terms 
in the variation of the bulk action. 

These boundary terms also allow one to derive an alternative and useful definition of the con- 
served currents. We do this by varying the off-shell action, now allowing a^(u — >■ °o) to vary, and 
then going on-shell by applying the equations of motion. Due to the boundary terms, only the 
surface term at u — > oo remains. Thus 



(34) 

on— shell 



In particular this relates the charge density in the boundary theory to the bulk radial electric field. 
The boundary term Sg(u m i n ) for ao should then be interpreted as a source term for this field. One 
could in principle also add boundary terms at u — > °°. These have no effect on the derivation of the 
equations of motion, but may change the value of the on-shell action, and therefore may lead to 



additional contributions to (34). 



A state with a non-zero baryon number density corresponds to an embedding with a radial elec- 
tric field. In particular, for a U embedding (in both the confining and deconfining backgrounds) this 
requires the addition of a baryonic source at the tip, corresponding to a uniform spatial distribution 



6 We would like to stress that this is a gauge invariant definition. The standard boundary condition on the gauge 
field in AdS/CFT fixes the value of cim{u —> °°). In this case only the transformations that vanish at u — > °o are gauged 
in the bulk. In particular, these transformations do not change the asymptotic value of oq. 
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of wrapped D4-branes (Figj6]l,e) [1481 (see also ||49l,l50l0 . We assume that the distribution is dilute 
enough so that we can ignore interactions between the D4-branes. We should therefore include the 
D4-brane action, which is given by (in the confining background) 

S D 4 = -jVVaudaNc Qw -oq (mo) J ) (35) 

where nm is the density of D4-branes. The first term is the D4-brane DBI action, and corresponds 
to the baryon mass, and the second term comes from the N c strings that connect each D4-brane 
to the D8-brane. The second term is precisely the boundary term at u m in = that was discussed 
above. 

The resulting asymptotic behavior of the gauge field is 

5 u 3 l L 

where d = N c no4 is the baryon number density. On the other hand, extremizing the action with 
respect to nm fixes the value of the gauge field at the tip to (wo) = uo/3 = m baryon /N c . This 
implies, as expected, that a non-zero density configuration exists only when the chemical potential 
is above the baryon mass. In fact the non-zero density state is always the dominant one. The 
transition to "nuclear matter" occurs at \Xy = mt, aryon /Nc. Near the critical point the density scales 
linearly with the chemical potential d ~ fly — m/, ar y on /N c - The D4-brane action also sources the 
embedding field x$(u), creating a cusp at u = uq. This can be understood in terms of a force balance 
condition between the D4-branes pulling down and the D8-brane pulling up. 

At high temperature the preferred embedding is parallel and there is a different finite density 
solution. In this case the gauge fields and d^ are independent, and the field theory has a con- 
served axial current as well as a baryon current. For a baryonic solution we take (°°) = jly 
and Oo(°°) = 0- m addition, since the D8-brane and anti-D8-brane reach the horizon we must im- 
pose <3q (ut) = % (ut) = 0. Therefore the radial vector electric field, and thus the baryon number 
density, is non-zero when fly > 0. In this phase d ~ T 3 /iv for small /ly. 



a b c d e 

Figure 6: D8-brane embeddings in the Sakai-Sugimoto model: (a) confined vaccum (b) deconfined 
vacuum (c) deconfined plasma (d) confined nuclear matter (e) deconfined nuclear matter. 
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3.3 Magnetic catalysis of chiral symmetry breaking 



A strong magnetic field in QCD is believed to catalyze the spontaneous breaking of chiral symme- 
try flU [2l [3J . The basic mechanism for this is that in a strong magnetic field all the quarks sit in the 
lowest Landau level, and the dynamics is effectively 1+1 dimensional. This phenomenon has been 
exhibited in the Sakai-Sugimoto model in [|5il l52|j^] 

With a uniform background magnetic field b, the D8-brane action in the deconfining back- 
ground becomes 



Sm=-<* J d 4 xduu 4 J (f(u){^{u)) z +^\ + • (37) 



The U embedding has the same form as before (26 ), but now uq depends on the magnetic field, as 



shown in Fig. [7^. The mass scale associated with chiral symmetry breaking is seen to increase with 
the magnetic field. One therefore expects that chiral symmetry breaking becomes more favored as 
the magnetic field increases. This is indeed the case, as can be seen by comparing the Euclidean 
actions of the U and parallel embeddings as the temperature and magnetic field are varied. The 
resulting phase diagram is shown in Fig.[7J). We observe that in this model the critical temperature 
approaches a finite value at infinite magnetic field. 

A qualitatively similar effect was observed in the D3-D7 model above. However in the D3-D7 
model there is a critical value of B/T 2 above which the chiral symmetry is always broken, whereas 
in the D4-D8 there is a critical temperature above which the chiral symmetry is always broken. 




Figure 7: Magnetic catalysis: (a) Chiral symmetry breaking mass scale, (b) Phase diagram. 



It is also instructive to study the effect of the background magnetic field on the mesons. This 
was partly done in [|55l . in which the high spin mesons were studied. It was shown that the mag- 
netic field enhances their stability by increasing their angular momentum, and thereby increasing 
the dissociation temperature at which they fall apart into their quark constituents. This is consistent 
with the above results. 



7 At non-zero baryon number density the magnetic field can actually induce an inverse magnetic catalysis in this 
model 15311541. 
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3.4 Anomalous currents 



The chiral anomaly leads to two interesting phenomena when both the magnetic field and chemical 
potential are non-zero. The first is the generation of anomalous currents in the chiral symmetric 
phase of QCD. In ll56l 1571 it was shown that the combination of a magnetic field and a non-zero 
baryon chemical potential generates an axial current 



Ja 



6 \lf ys B. 



2n 2 



(38) 



Since the source for this current is the anomaly it is an exact result, and should be valid in particular 
at strong coupling. Similarly, an anomalous vector current is generated in a non-zero axial chemical 
potential: 

e 



Jv 



phys 



B. 



(39) 



This is known as the "chiral magnetic effect", and may have some relevance to heavy ion physics 
at RHIC [5J. Within the D3/D7 model, this effect has been discussed in 11551 . 

In the Sakai-Sugimoto model, the chiral- symmetric phase corresponds to the parallel D8- 
D8 embedding in the deconfined background, and the chiral anomaly is encoded in the five- 



dimensional CS term (24). The background magnetic field and chemical potentials correspond 
to different components of the worldvolume gauge field. Through the five-dimensional CS term 
these source a third component, which corresponds to a current in the four-dimensional theory 
[59|. Let us review the calculation of the anomalous axial current in the Sakai-Sugimoto model. 
The calculation of the anomalous vector current is virtually identical. 

To be specific, we will consider a background magnetic field in the x\ direction by turning on 
a background gauge field = X2b. A non-trivial boundary value of will then source, via the 
CS term, a non-trivial of. In general and of can depend on both u and X2 in this case, although 
on-shell they will depend only on u. We will take <3q (°°) = an d «f (°°) = 0. The D8-brane DBI 
and CS actions in this case become 



S d D % = £d 4 xduu 5 ^ 2 J(l-(al') 2 + f(u)(c4')^ + (40) 



Scs = — / d 4 xdu b{a^a^ — a 'af ) + (^(a^' d2d\ _ e^o^i 



The necessary boundary terms are 



*2- 



*2~ 



(41) 



(42) 



By integrating by parts the last two terms in the CS action one can show that up to a surface term 
at u — Y °° , the bulk CS and boundary actions combine into a bulk action 



Scs + $d = J d 4 xdu -b(a,Qa^' — a^'af ) — -aj \a,Q t^af — c^ag af ) 



(43) 
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One can get rid of the remaining surface term by adding a boundary term at u — > °° 115911 , however 
this particular term does not contribute to the on-shell action, so we might as well ignore it j^] 
The equations of motion for Oq{u) and a\ (u) can be integrated once to yield 

y/u 5 + b 2 u 2 aY,' (u) A , , 

- -3bcrf(u)+d (44) 



i-( a v'( u ))i+f(u)(c4'(u)y 

= = -3bal(u), (45) 



Vu 5 +b 2 u 2 f(u)af(u) 



^\-(a v \u)) 2 + f(u)(af(u)) 2 
where d is the baryon number charge density. The integration constant in the a\ equation vanishes 



since {uj) = and f(uj) = 0. Using (34) and (45) we can then evaluate the axial current: 



j l A = lba v H = ^ v . (46) 

The correctly normalized physical currents are given by J = l^llta' ,yV /R 5 )j, where the factor 
of 2 comes from adding the anti-D8-brane contribution, and the physical chemical potentials are 
jjphys _ (r/(2tiu'))ii. Thus in terms of the physical variables our result translates to 

T A = ^rf»B. (47) 
Similarly, for an axial chemical potential we would get 

Jv = ^ hys B. (48) 

Interestingly, our results are half of the weak-coupling results (where e = N c in the holographic 
model). It has been argued that the discrepency in the axial current is due to a different treatment 
of the triangle anomaly, consistent vs. covariant, and can be corrected by adding an appropriate 
Bardeen counterterm on the boundary lloTll . However, the same counterterm leads to a vanishing 
vector current. This issue is still under investigation. 

3.5 The pion gradient phase 

In the broken chiral symmetry phase the chiral anomaly leads to a novel finite density phase that 
dominates over nuclear matter at large magnetic fields [4]. In this phase the baryon charge is 
carried not by baryons but rather by a non-zero pion gradient background: 

D =4^^°' <49) 



Other boundary terms at u — >• °° could affect the on-shell action, and therefore the currents. See for example 
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where 



i? hys B. 

In the Sakai-Sugimoto model (in the a u = gauge) the pseudoscalar meson appears in the zero 



w (50) 



mode of (27), so 



d^(p(x^) = -a^(x^,u^oo). (51) 

As in the chiral- symmetric phase, the presence of a vector chemical potential together with a back- 
ground magnetic field sources a component of the axial gauge field, which in this case corresponds 
to a non-trivial gradient of the pseudoscalar field [|59ll62l . Since there is only one flavor this field 
should really be thought as the r\' meson. For simplicity, we will consider only the confined phase 
with the anti-podal D8-brane embedding, namely uq = ukk- (The results do not change qualita- 
tively for more general U-shape embeddings, or for U-shape embeddings in the deconfined phase.) 
Following 11591 , the D8-brane DBI action in this case is 



s^ /= -^^2^ m " 5/ v(^-w 2 +(<) 2 ) (52) 



and the CS plus boundary actions are the same as in the deconfined phase (43 ). The equations of 
motion integrate to 



y/u 5 + b 2 u 2 Oq'{u) 



^-(a >))2 + (a>))2 

Vu 5 +b 2 u 2 af(u) 
v / 7 ^-(a ^W)2 + (af( M ))2 



-3bc4(u)+N c n D4 (53) 



-3teo(«)+c, (54) 



where we have explicitly included the baryon sources in the cIq equation. Our boundary conditions 
are now Oq (°°) = /iy and a\ (°°) = —V(p(x^). In particular a\ (°°) is a field rather than a parameter 
in the boundary theory, and we must minimize the action with respect to its value. This simply sets 
j\ = and therefore sets the integration constant in the of equation to c = ^pily^ 

We can now compute the total baryon number charge density d using the same procedure as in 
the previous section for the current. In the absence of sources hda = and we find 

d = -^bc4H = 3 -bVcp. (55) 
Let us express this in terms of the physical variables. First we must define a field with a canonically 



normalized kinetic term. Inserting (27 ) into the action (52 ) we find that the canonically normalized 
field is given by 



This is also consistent with the fact that there are no quarks in this phase to carry such a current. 



19 



Converting to physical variables we then find 



D= A^ S ^'< <5?) 



in agreement with (49). We would like to stress that this agreement did not depend on the specific 
value of f„i required for canonical normalization, since it cancels out when we express the result 
in terms of (p. The correct numerical factor of l/(4n 2 ) is a direct consequence of including the 
proper boundary terms in the action, leading to the "3/2" in p5|). 



To find the value of the gradient V<p we need to solve (53) and (54). The result will not be 
as simple as (50). In particular it is not linear in the magnetic field, since we are using the full 
non-linear DBI action. It turns out that a closed form solution can be found in terms of a new 
variable 

f" 3bdu 

y= / =■ (58) 

JllKK 

The solution is 



yfmJuS (1 +*2fi-3) - (IbHv) + (3W<p)2 



v nv fcoshy \ A n v sinhv 

where y x = y(u — > °°). The pseudoscalar gradient is then given by 

V<p = -of (oo) = ^ tonhy„ . (60) 

The dependence on b is shown in Fig. [8^. For small b the behavior is linear in b: 

V<p« -^n v b, (61) 
2u KK 



and in terms of the physical quantities: 



in agreement with the single flavor version of ( |50| ). 

As in the case with no magnetic field, an embedding that includes sources is possible above 
a critical value of the chemical potential, which then becomes the dominant configuration. This 
describes a "mixed phase" that includes both "pion-gradient" matter and nuclear matter, with a 
total baryon number density 

d=^bV(p + N c n D4 . (63) 

As before, one can find a closed form solution for and of, and from it determine the values 
of V<p and no4 in terms of the magnetic field b and baryon number chemical potential /ly. The 
resulting phase diagram is shown in Fig. [8J5. In particular, the critical value of /ly is determined 
by setting = 0. The relative proportion of baryons in the mixed phase increases with fly and 
decreases with b. 
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Figure 8: a. The pion gradient, b. Phase diagram with magnetic field and baryon chemical 
potential. 



3.6 Magnetic phase transition 

The high-temperature chiral-symmetric phase of this model also exhibits an interesting magnetic 
phenomenon associated with the distribution of the baryonic charge [63J. In what follows we ana- 
lyze the situation for the D8-brane. As in the broken chiral symmetry phase above, the distribution 
of baryonic charge along the radial direction u changes with b. We can therefore identify two types 
of baryonic charge, one originating from the horizon, and the other from outside the horizon. The 
latter d*, corresponds to D4-branes that are radially smeared inside the D8-brane. This can be best 
seen from the longitudinal and transverse conductivities ll64l 



\ u\ -\-b 2 Uj + wL(d — d*) 2 

°L = 1 (64) 

u^ + b 2 

b(d — d*) d* 

ul + b 2 b 



Ot = v 3 , ,;' +-^, (65) 



where d* = —3ba\ (uj) ■ In particular only the horizon charge d — d* contributes to the longitudinal 
conductivity. In the transverse conductivity, the horizon charge contributes as an ordinary dissi- 
pative fluid, whereas the charge outside the horizon d* behaves as a dissipation-free fluid. This is 
consistent with an interpretation of d* as the charge filling the lowest Landau level. As the mag- 
netic field increases more of the charge is "lifted" from the horizon, representing the transition to 
the lowest Landau level in the boundary theory. 

In fact for a fixed density at low enough temperature this transition is a first order phase tran- 
sition as a function of the magnetic field, in which the charges jump into the lowest Landau level. 
This is easiest to see in the zero temperature limitf^] In this case one can solve the gauge field 
equations analytically in terms of a variable 

3bM (66) 



o Ju 5 + b 2 u 2 + d 2 cosh 2 z, 



10 Strictly speaking, at zero temperature the theory is in the confining (and broken chiral symmetry) phase. We are 
considering the meta-stable state obtained by adiabatically reducing the temperature. 
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where = z(u — > °°). The solution is 

v dsinhz A Jcoshz d 
B^coshzoo ' 1 3^coshzoo 3b ' 

There are actually three solutions, representing two stable phases and an unstable phase. As the 
magnetic field b is increased, for a fixed total baryon number density d, one finds a first order phase 
transition between the two stable phases. Both the magnetization and the chemical potential are 
discontinuous in this transition (Fig. [9]), which is reminiscent of a metamagnetic phase transition. 
The large magnetic field phase represents the situation where all the charge is in the lowest Landau 
level, with the chemical potential and free energy given by 




Figure 9: (a) /i and (b) M as functions of b for d = 1 and T = 0.09, below the critical point. 
There are now two branches of stable solutions, and the phase transition between them occurs at 
£ = 0.235. 



The magnetic transition persists also at non-zero temperatures that are low relative to the den- 



sity d. Too a very good approximation this happens when b ~ d 2 ^ 3 (Fig. 10 1), which is the behavior 
expected for the lowest Landau level. At high temperature the transition disappears (Fig.[l0p). 





Figure 10: (a) The phase diagram in the d-b plane at T = 0.07, and (b) the critical line in the T-d 
plane. 
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4 The D3-D7' model 



The study of magnetic properties of planar matter is a very active area of research in condensed 
matter physics. A simple holographic model for charged fermions in three dimensions can be 
obtained by T-dualizing the D4-D8 setup. This leads to a D3-D7 configuration with the two sets 
of branes intersecting on a plane lfT2l [T4l (for a related model using a D2-D8 system see [65 1). 
However, unlike in the D4-D8 configuration, here the branes have a mutually transverse coordinate. 
On the one hand this allows the fermions to be massive, but on the other hand it leads to an 
instability since the different branes repel. 



4.1 Stable embeddings 

First we have to address the issue of stability. As before, we will employ the probe approximation 
and consider a single probe D7-brane. The background (at finite temperature) in this case is 

L- 2 ds\ Q = r 2 (-h(r)dt 2 + dx 2 + dy 2 + dz 2 )+r- 2 ( ^— + r 2 dQ.i] (69) 

\h(r) J 

F 5 = 4L 4 (r 3 dt A dx A dy A dz A dr + dD. 5 ) , (70) 

where h(r) = 1 — r 4 /r 4 , rj = nLT and L 2 = y/4ng s N c (x'. It is convenient to parameterize the 
five- sphere as an S 2 x S 2 fibered over an interval: 

dQ. 2 5 = dy 2 + cos 2 ^{d^f + sin 2 ^(dQ.^) 2 , (71) 

where < i/a < n/2. The first S 2 shrinks at the "south pole" iff = Tt/2 and the second S 2 at 
the "north pole" \y = 0. The D7-brane wraps the two S 2 's and extends along (x,y), and has an 
embedding described by z{u) and i//(m). In particular y is dual to the fermion bi-linear operator in 
the field theory corresponding to the fermion mass. The D7-brane DBI action in this background 
is given by 

S DB1 = -4,yy J d 3 xdrr 2 cos 2 ^sin 2 y/^ 1 + r 4 h(r)z' 2 + r 2 h(r) \\f' 2 , (72) 

where ,jY = 47T 2 /l7L 8 / 'g s . A massless embedding would correspond to ty/ = n/A. However the 
fluctuations contain a mode that violates the Breitenlohner-Freedman bound, and therefore the 
embedding is unstable. This can also be seen by trying a more general embedding with a large r 
behavior of 

V/(r)~| + cr A . (73) 

The equation of motion for i// gives A(A + 3) = —8, which does not have a real solution. 

Fortunately, the D7-brane can be stabilized by turning on some worldvolume flux 11661 , in this 
case on the two- spheres [fT51 : 

2%o!f = y {f x ddp + / 2 ^ 2) ) , fi= l7 ^f- m (m e z) . (74) 



23 



This changes the DBI action, 



Sdbi 



/ d 3 xdrr 2 J (4cos 4 V/ + /i 2 ) (4sinV + / 2 ) {I + r*hz' 2 + r 2 h\j/' 2 ), (75) 



and there is now also a CS term which gives, 




(76) 



The asymptotic behavior of \jf(r) is now 



(77) 



where is determined by the solution of 



(fi + 4 cos 4 y/oo) sin 2 \ffoo = (/| + 4 sin 4 y/oo) cos 2 y/ c 



(78) 



and 



/ 2 + 16 cos 6 y/oo - 12 cos 4 y/ c 



/ 2 + 4cos 6 v/oo 



'oo 



(79) 



In particular, the embedding is stable for a large enough flux. The coefficient of the leading term is 
related to the fermion mass, and that of the subleading term corresponds to the bi-linear condensate. 
Note that the scaling dimension of the bi-linear operator is given by — A . This represents a large 
anomalous dimension, which is not surprising given that the model is non-supersymmetric. We 
should require however that the operator be relevant, namely that A > —3, and therefore that 
A+ < 0, in order to consider the leading term as a "mass deformation". 

Generally there are two types of embbedings, that differ in their small r behavior: Minkowski- 
like (MN) embeddings, in which the D7-brane terminates smoothly outside the horizon (Fig.fTTfc), 
and black-hole (BH) embeddings, in which the D7-brane crosses the horizon (Fig. [Tip). We refer 
the reader to lfT5l for the explicit embedding equations for \j/(r) and z(r), and for their numerical 
solutions. 

In an MN embedding y(ro) = 71 / 2 or for some tq > rj, corresponding to one or the other S 2 
shrinking. This indicates that the dual field theory has a mass-gap related to r§ — rj. An important 
condition for the existence of MN embeddings is the absence of sources for the worldvolume gauge 
field. Unlike in the model of the previous section, there are no localized sources in this model. All 
sources correspond to branes (or strings) connecting the D7-brane to the horizon. These inevitably 
pull the D7-brane down to the horizon, resulting in a BH embedding instead. This means that the 
flux on the S 2 that shrinks must vanish. For yA(ro) = 7C/2, which means that f\ = 0, we find stable 
massive embeddings in the range 



0.5235 < y-, < 0.6251 . 



(80) 
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(There are also the "mirror" embeddings with y(ro) = and f% = 0.) Note that the allowed 
values of ty/oo are quantized since the stabilizing flux is quantized. There are also two isolated MN 
embeddings with = and y^o = Jc/2. 

BH embeddings describe gapless phases in the dual theory. These embeddings exist generically 
for any f\ , satisfying the stability condition. 



r sin ip r sin V 




r cos r T r cos i/> 



Figure 11: (a) An MN embedding with y(ro) = 7t/2. (b) A BH embedding. 



4.2 Finite density and background fields 

For embeddings corresponding to finite density states in a background magnetic field we need to 
turn on the appropriate components of the worldvolume gauge field. As in the D4-D8 model we 
will work with the dimensionless field = (27ta' /L)A„. There are additional terms in the DBI 
action, 

Sdbi = j drr 2 ^ (4cos 4 V + ~ff) (4sin 4 y + / 2 2 ) x 

X J(l + r*h{r)z' 2 + r 2 h(r) y' 2 - a' 2 ) ^1 + ^ , (81) 

and also in the CS action, 

Scs = -JTfxh J drr 4 z'(r) + 2yK J drc(r)ba' (r) , (82) 

where 

If 1 1 

cW = ^Z? / 5 2 x5 2 C4 ^ r)) = V^)--sm4i/<r)-i//oo + -sin4i^ (83) 

We have fixed a gauge for the RR field such that c(°°) = 0. For MN embeddings we must also add 
a boundary term at r = tq (as explained in section pjrj 



S d (r )=2<yKc(r )bao(r ). (84) 



In ifTSIl this term was derived by demanding invariance of the CS term / C4AFAF under gauge transformations 
of the RR field and then fixing c(°°) = 0. However it can also be obtained by canceling the surface term in the variation 
of the CS term, when we present it as J F5 A A AF. 
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The quantity c(ro) has a nice physical interpretation: it is the total amount of 5-form flux captured 
by the D7-brane in the MN embedding. It is completely fixed by the asymptotic value of the 
embedding angle For BH embeddings the boundary term vanishes since ao(/r) = 0. 
The integrated equation of motion for ao(r) is given by 



G(r)a' (r) = d-2bc(r), 



(85) 



where 



G{r) 
Y{r) 




(/ 1 2 + 4cos»(/ 2 2 + 4sin» 
Y(r) 



+ {l+hr 4 z' 2 + hr 2 xi/ 2 -a 2 ) 



(86) 
(87) 



and d is the total charge density. As in section[3j we are using (34) to define the conserved currents. 
The quantity on the RHS of (85 ), d(r) =d — 2bc(r), is the contribution to the charge density from 
radial positions below r. 

We would also like to study the response of the system to a background electric field. To 
this end we should consider a more general ansatz for the gauge field with a x (t,r) = te + a x (r), 
a y (x,r) — xb + a y (r), in addition to ao(r). The current densities will be contained in the asymptotic 
behaviors of a x (r) and a y {r). The gauge field equations in this case become 



G(r) a' (r) 
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(88) 

(89) 
(90) 



where j x is the longitudinal current density and j y (r) is defined by analogy with d(r) as j y (r) 



j y — 2c(r)e, where j y is the transverse current density. The factor G(r) is defined as before (86), 
now more generally with 



-( 1 + 74)«0+^+(l-^4)^. 



2eb , , 

~r a Q a y 



(91) 



4.3 Quantum Hall states 



Let us consider first the response of the gapped MN embeddings. This is determined by the re- 
quirement that there are no sources, namely by regularity of the gauge field at r = tq. This implies, 
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in particular, that 



d(r )=d-2c(r )b = 0. (92) 

The entire charge in the MN embedding is thus due to the CS term and corresponds to a "fluid" of 
D5-branes inside the D7-brane. The charge density is proportional to the magnetic field, and the 
proportionality constant is fixed by the value of c(ro), and therefore of This is the key property 
of a quantum Hall state, which is characterized by a specific quantized value of the Landau-level 
filling fraction v d/b. In terms of the physical variables D = (2na'^V /L 4 ) d and B = b/ (2na'), 
the filling fraction is given by 

2%D 2N C , . 
V =— = -±c(r ). (93) 

For the range of values of i/4o needed for stability ([80]) we get 



0.6972 <^< 0.8045. (94) 
Furthermore, the filling fractions are quantized according to the quantization of 1//^. The actual 



numbers can be obtained by solving (78), for a specific flux f2 (with f\ = 0), and plugging into 
( 83 ) with vt^o) = ^/2, but they are not particularly illuminating (for example, they are not rational 
numbers). The isolated embeddings with i/4o = and 1//^ = n/2 correspond to v/N c = 1 and 0, 
respectively. 

The current densities can likewise be computed by requiring regularity of the spatial compo- 
nents of the gauge field. This condition implies that 

j x = and j y (r ) = j y - 2c(r )e = , (95) 

from which we can deduce the longitudinal and transverse conductivities: 

v 

2% 



0, C7xy=— . (96) 



Thus the MN embeddings, when they exist, describe quantum Hall states with quantized transverse 
conductivities, and vanishing longitudinal conductivities. Furthermore, in the holographic descrip- 
tion, the quantization is topological since it originates from the Dirac quantization of the magnetic 
fluxes on the S 2 's. In particular, Oxy in the MN embeddings is independent of the temperature. 

Quantum Hall states are gapped to both charged and neutral excitations. In this model charged 
excitations are described by strings stretched from the D7-brane to the horizon, and therefore have 
a mass proportional to r$ — rj. This is seen to increase with the magnetic field, as shown in 



Fig. 12 1. The neutral excitations correspond to fluctuations of the D7-brane worldvolume fields, 
and are also found to be massive [67J (see also ll68ll ). The spectrum of neutral excitations includes 
a magneto-roton, which is a collective excitation whose dispersion relation has a minimum at non- 



zero momentum (Fig. 12 i). A similar phenomenon is seen in real quantum Hall states H691 
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Figure 12: (a) Mass-gap for charged states as a function of the magnetic field, (b) Dispersion 
relation of the two lowest neutral modes, showing the magneto-roton. 



4.4 Fermi-like liquid 

BH embeddings describe gapless Fermi-like liquids. For a BH embedding, d{rj) corresponds to 
the horizon charge density carried by the "quarks", and it need not vanish. In particular, if we add 
sources to an MN embedding, thereby violating ( [92] ), it deforms continuously into a BH embedding 
with horizon charge. 

To compute the electrical response in a BH embedding we have a couple of options. The 
standard approach is to extract the conductivities using linear response from the current-current 
correlators, computed holographically by studying fluctuations of the bulk gauge fields to quadratic 
order. The other option is to find a consistent solution in the presence of an external electric field 
|70l. The advantage of the second approach, when it is applicable, is that it gives the complete 
non-linear response. Using this method for the BH embeddings one finds 



Nr. 



lit 1 \b 2 + r 



-d(r T ) + 2c(r T )) (97) 



T / 

° xx = ^2 ^4 yfiM 2 + iff + 4cos4 YM) {fl + 4 sin 4 w(r T )){b 2 + 4). (98) 

Note that the transverse conductivity has two components. The first involves the horizon charge, 
and resembles the contribution of an ordinary dissipative system of charges. The remaining charge, 
corresponding to the fluid of D5-branes inside the D7-brane, contributes like a dissipationless 
system. The longitudinal conductivity involves only the first component. This is basically the 



same separation that was seen in the D4-D8 model (see equations (64), (65)). 

This state of holographic matter exhibits a variety of other interesting phenomena as a function 
of the charge density, temperature, background magnetic field, and mass. 

Consider first the state at T = 0, d = 0, b = and m = 0. In this case the D7-brane embedding 
is actually AdS$ x S 2 x S 2 , so this situation is described by a conformal field theory. Note that in 
this case o xx ^ 0. At non-zero density the system becomes unstable to the formation of stripes. The 
instability is signaled by the existence of a quasi-normal mode (in the transverse gauge field sector) 
with positive imaginary part in a finite range of momenta 11711 . At a high enough temperature or 
high enough magnetic field this instability disappears. It is convienient to parametrize the situation 
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with b = b/r^ and d = d/r^. Then the instability towards the striped phase at zero magnetic field 
and m = happens for d > 5.5. This is demonstrated in Figure 13 1. Above the quantum critical 



5 r 



T 




10 20 30 40 50 



Figure 13: (Left)- Phase diagram in the T-d plane showing the quantum critical point. (Right)- The 
phase diagram in the b, d plane for m = 0. Above the line the system is in the homogeneous phase 
and below the line in the striped phase. 

point (T = m = d = 0) there is a region which resembles a Fermi-like liquid, and on both sides 
there is a striped phase. For non-zero b and m = the instability sets in at some other value of d, 
as shown in Fig. [T3^ > [72] . As m increases the instability sets in at a lower temperature. 

The system also has a zero sound mode. At non-zero temperature the quasi-normal mode 
with the smallest imaginary part at low momentum is a purely dissipative hydrodynamical mode 
(co(k = 0) = 0). At some non-zero momentum it meets another purely dissipative mode (ft) purely 
imaginary) and crosses from a hydrodynamical regime into a collisionless regime, where the re- 
sulting complex mode can be identified with the finite temperature zero-sound mode [71] (see also 
11731 ). The zero sound mode becomes massive as the magnetic field crosses a critical value [172117411 
(FigfT^i). 



0.8 „ M 




Figure 14: (Left)-The line in the b,d plane for m = 0, separating the situation when the zero- 
sound is gapped (above the line) and the situation when it is gapless (below the line). (Right)-The 
magnetization as a function of magnetic field (solid line) and the individual contributions from the 
DBI (dashed line) and the CS term (pointed line) for m = 1 d = 1 and rj = 0.1. 
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For »i^0 the system can have a non-zero transverse conductivity, even at zero magnetic field. 
This is due to having a non-trivial c(r) for these embeddings. This is reminiscent of the anomalous 
Hall effect (AHE) that appears in ferromagnetic materials (for a review see [?]). Indeed for m ^ 
the system is ferromagnetic (Fig. 14?) due to the second term in (82). Note that both the AHE 
and the ferromagnetic behavior, as well as the instability towards a striped phase, have a common 
origin in the Chern-Simon term J drc(r)F AF in the brane action. 
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